Extended Superscaling of Electron Scattering 

from Nuclei 



C. MaieronQ, T. W. Donnelly 
Center for Theoretical Physics, Laboratory for Nuclear Science 
and Department of Physics 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139-4307, USA 

and 

Ingo Sick 

Departement fiir Physik und Astronomic, Universitat Basel 
CH4056 Basel, Switzerland 



An extended study of scaling of the first and second kinds for inclusive 
electron scattering from nuclei is presented. Emphasis is placed on the 
transverse response in the kinematic region lying above the quasielastic 
peak. In particular, for the region in which electroproduction of reso- 
nances is expected to be important, approximate scaling of the second 
kind is observed and the modest breaking of it is shown probably to 
be due to the role played by an inelastic version of the usual scaling 
variable. 

February 8, 2008 



1 Introduction 

In recent studies the concepts of scaling of the first and second kinds and 
superscaling have been explored, focusing on the region of energy loss at or below 
the quasielastic peak in inclusive electron scattering from nuclei. Scaling of the 
first kind corresponds to the following behavior: if the inclusive cross section is 
divided by the relevant single-nucleon electromagnetic cross section {i.e., weighted 
by the proton and neutron numbers, Z and N, respectively, and with appropriate 
relativistic effects included — see for detailed discussions), then at sufficiently 
high values of the momentum transfer q the result so-obtained becomes a function 
of a single scaling variable and not independently of q and the energy transfer u. 

-'^ Present address: Fachbereich Physik, Universitat Rostock, D-18051 Rostock, Germany 
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Various definitions for tlie scaling variable exist (see, for instance, [Q); however, 
when q is high enough they are almost always simply functionally related in 
ways that yield scaling behavior in all cases. Such g-independence of the reduced 
response F{q,uj) —* F{iIj), where ip = ip{q,uj) is the scaling variable, is called 
scaling of the first kind. Moreover, motivated by earlier work [Q, it has been 
found that when the typical momentum scale of a given nucleus k is appropriately 
incorporated in the definition of the scaling variable and the reduced response is 
also appropriately scaled, F{q,u) f{<i,^) = kx F{q,Lj), then a second type 
of scaling behavior is seen — the result becomes independent of nuclear species. 
This is called scaling of the second kind. When scaling of both the first and 
second kinds occurs, one calls the phenomenon superscaling. 

The studies undertaken recently |^, 0] demonstrated the quality of the scaling 
behavior, finding that scaling of the first kind is reasonably good for ip < (below 
the quasielastic peak), that scaling of the second kind is excellent in this region, 
but that both are violated for > (above the quasielastic peak). Indeed, the 
scaling of the first kind has been known for some time to be badly violated in this 
latter region. The recent analyses showed that so also is scaling of the second kind 
broken in the ijj > region, but much less so. Finally, in p| an initial attempt was 
made to use the limited information on the separate longitudinal and transverse 
responses, and hence on the scaling behavior of the respective reduced responses. 
The former appears to superscale, whereas the dominant scale-breaking effects 
appear to reside more in the latter. 

In the present work we pick up these ideas and extend them. We begin by 
updating our analysis of the relevant inclusive electron scattering scaling behavior 
in the ip < region, thereby obtaining refined values for the typical nuclear 
momentum scale (henceforth, as in our previous work which was motivated by 
the Relativistic Fermi Gas (RFG), called the Fermi momentum kp) and for a 
small energy shift Eghift included to have the quasielastic peak occur at the place 
where the scaling variable is zero. Moreover, we assess the sensitivity of the 
results to variations in both kp and Eghift to provide some idea of how much 
change in one or both can be tolerated when studying the region where ip > 0. 

We wish to focus on the scale-breaking effects, especially in the ip > region, 
and accordingly we have isolated the transverse response using our previous ap- 
proach. Within uncertainties that are unfortunately not as small as is desirable, 
and only for a limited range of kinematics, it appears that the longitudinal re- 
duced response in fact superscales and reasonably satisfies the Coulomb sum 
rule. We shall assume that this is a universal behavior (for all kinematics and 
for all nuclei — that is the impact of having superscaling), shall then remove the 
longitudinal contributions and thereby obtain the transverse reduced response, 
fx- Naturally the uncertainties in knowledge of fi propagate into corresponding 
uncertainites in /y, and given better L/T separations the procedure could be 
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refined. However, it is important to the best of our current ability to isolate the 
transverse part of the inclusive response as it is the one expected to contain the 
leading scale-breaking effects 0. 

Once fx has been isolated we focus on the region above the quasielastic peak, 
exploring the scaling behavior of Jt as a function of q and kp, namely, for first- 
and second- kind scaling. We further divide our discussions into two regimes, (1) 
a study of the resonance region where A's and A^*'s are expected to play an 
important role along with non-resonant meson production, and (2) a study of the 
very limited data available in the Deep Inelastic Scattering (DIS) region. Upon 
seeing that the second-kind scaling behavior is only moderately violated in these 
regimes, and motivated by the type of analysis performed in studying the EMC 
effect, we also define an appropriate ratio involving the /j-'s for a pair of nuclei. 
This provides a convenient measure of the extent to which scaling of the second 
kind is or is not respected. Indeed, as we shall see, the ratio is very close to 
unity from the most negative values of the scaling variable {i.e., far below the 
quasielastic peak) up through the resonance region. Only in the DIS region does 
the ratio differ from unity by as much as 25%; for most of the region the results 
lie typically within about 10% of unity. 

The paper is organized as follows: in the next section we briefly summarize 
the basic scaling formalism, drawing on our previous work JI], 0] where more 
detailed discussions can be found. In Sec. |] the updated determinations of kp 
and E shift are discussed, in Sec. ^the transverse scaling functions fx are presented 
and in Sec. H the ratios involving pairs of nuclei are introduced. For the last, the 



discussions are focused on two kinematic regimes, the resonance region in Sec. 5.1 



and the DIS region in Sec. |5]^. In the former additional modeling is presented to 
help in understanding why the ratios differ from unity even by the small amount 
they do. Finally, in Sec. ^ we summarize our observations and conclusions. 



2 Basic Scaling Formalism 

We begin by summarizing some of the essential expressions used in previous 
studies of scaling and superscaling in the quasielastic region. First, using the 
nucleon mass ttiat as a scale, it proves useful to introduce dimensionless variables 
to replace the 3-momentum transfer q and energy transfer tu, namely 

A - ^ (1) 

2mAr 

The (unpolarized) inclusive electron scattering cross section then depends only on 
K = a and the electron scattering angle 9^- The dimensionless 4-momentum 
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transfer squared is given by 



T=\Q'\lAm'^ = K'~\\ (3) 

and = cu^ — < in the conventions used here. 

In past work it became apparent that a convenient dimensionless scale in 
quasielastic electron scattering from nuclei is provided by the ratio of a charac- 
teristic nuclear momentum k = \J <k'^ > to the nucleon mass m^. For example, 
in the Relativistic Fermi Gas (RFG) model the characteristic momentum is the 
Fermi momentum hp and the dimensionless scale is given by r/i?, where 

^ « 1, (4) 

where typically the Fermi momenta range from as small as 55 MeV/c for deu- 
terium, 200 MeV/c for ^He, to as large as about 250 MeV/c for very heavy 
nuclei, and as a consequence the strong inequality above holds. A corresponding 
dimensionless energy scale is also useful. 

Bp = ^l + rfp = 1 + + . . . , (5) 

where then 

iF = ep-l = ]^'nl + ---. (6) 

Naturally the RFG is only a first approximation to the nuclear dynamics involved 
in the quasielastic region and thus the Fermi momenta actually employed — and 
these are obtained by fitting data as discussed below and in ^ — should 
really be regarded as effective parameters in the problem. Clearly across the 
periodic table the densities of nuclei change and this is reflected in the fact that 
the characteristic momentum scale (here kp) should also vary, roughly such that 
the density is proportional to k'p. As a consequence the width of the quasielastic 
response goes as kp. 

In past studies of the region at and below the quasielastic peak it has proven 
to be very useful to introduce scaling variables (see, for example, 0). In the 
most familiar approach the y-scaling variable is employed and one finds that at 
high momentum transfers the experimental inclusive cross sections divided by an 
appropriate single-nucleon cross section scale, i.e., become functions only of y 
and not independently of the energy or momentum transfer. Such behavior is 
called scaling of the first kind. Alternatively, again using the RFG for guidance, 
a dimensionless scaling variable ip emerges naturally 

■1' - , (7) 
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At the naive quasielastic peak where A = r (which corresponds to u = \Q^\/2mN) 
one has ip = and finds that the RFG response region is mapped into the range 
— 1<4'<+1. While appearing at first sight to be quite different, in fact the 
variables y and ip are closely related. One can write M 



y{Es = 0) 



l + 0{vF^,y/Ml 



(8) 



where Eg is the separation energy, the difference between the sum of the nucleon 
plus ground-state daughter masses and the target ground-state mass. Up to the 
choice made in making the scaling variable dimensionless by dividing hj kp (see 
below), when Eg is set to zero the two variables differ only at order rjp and order 
[M^_^]~^, and accordingly when scaling of the first kind occurs in one it is bound 
to occur in the other, as long as corrections to the leading order expressions are 
small. In effect, what the conventional y-scaling variable does that ip does not 
is to take into account the (small) shift in energy embodied in Eg. The simple 
relationship begins to fail when large excursions are made away from the QE 
condition ip = — in the present work the full expressions are always used. 

From such considerations one sees that an improved phenomenological di- 
mensionless scaling variable can be employed in treatments of superscaling [Q, 
namely one with an empirical shift E shift- This was done in previous analyses 

introducing a dimensionless scaling variable as above. 



X' 



n + x' 



r + K\ T [r 



(9) 



where XsUft = EsMftl^niN, A' = A - Xshift and t' = - X'' 

After this brief summary of the conventional choices of scaling variables, let us 
next turn to the scaling functions. We begin with the inclusive electron scattering 
cross section itself, which may be written in various forms (see also later): 



d^a / dVL(,duj 



+ St 

(Jm [vlRl{i^, X) + vtRt{i^, A)] 



W2 + 2Wi tan^ 9e/2) 



(10) 



where the familiar electron kinematical factors in this Rosenbluth form are given 
by 



Vt 



T 

V. 

r 



2k' 



+ tan^ 9^/2 



fill 
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and the longitudinal (L) and transverse (T) response functions are related to Wi^2 
via 



Rt 
Rl 



2Wi 

212 



K 
T 



W2 - —Wi. 

T 



;i2) 



The strategy in discussing scaling in the quasielastic region is, to the extent 
that it is possible, to divide out the single-nucleon eN elastic cross leaving only 
nuclear functions. As discussed in |, this can be accomplished using the 
reduced response 



(7m [vlGl{i^, A) + vtGt{k,, A)] 



(13) 



for the total cross section, or, when individual L and T contributions are being 
considered (as they are in part of the present work), using 



Fl = 
Fx = 

Here the functions Gl^t are given by 

Gl(«,A) = 

GTiK,\) = 

which involve the function A: 



Rj 



Gl{k, A) 
Rt 



Gt{k, A)" 



2«:[l+^^(l+^2)/2] 
2 

M 



2tGI, + 



2/c[l + e^(l + ^2)/2] 



A 



eF(i 
1 



r(l + r 



K 



1 - + o[v'f]. 



As usual one has 





= ZG%p 




= Z G\ip 


Wi 


= tGIi 




1 


W2 






1 + r 



■y2 

'En 
-y2 



Ge 



M 



(14) 
(15) 

(16) 
(17) 



(19) 
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involving the proton and neutron Sachs form factors GEp,n and GMp,n weighted 
by the proton and neutron numbers Z and A^, respectively. 

In the region of the QE peak where is small it is a good approximation to 
set A to zero and to expand the functions Gl,t in powers of r]p, retaining only 
the lowest-order terms, namely to use 

GL{n,\) = ^Gl + 0[vl] (20) 
Zt 

Gt{k,X) = -Gl + 0[7]l]. (21) 

K 

While such leading-order expansions in rjp are very good when r/i? x |?/^| is small, 
they become less so when very large excursions away from the QE peak are made 
and accordingly in the present work we have always used the full expressions in 
Eqs. dlMl- 

As discussed in 0, ^ the RFG model then yields scaling of the first kind, 
namely, the F's become functions only of ip (independent of k, that is, of q); 
indeed, as discussed in [|I], 0, so do the data when ip < 0. Moreover, as also 
discussed in ||, the RFG model and the data both display scaling of the 
second kind in that the F's can be made independent of the momentum scale in 
the problem — that is, independent of kp to the order considered in the expansion 
in the small dimensionless parameter riF- This is accomplished by defining 

/ = kpx F (22) 
U^T = kpxFL,T, (23) 

namely by making them dimensionless through multiplication by the factor kp. 
Comparing Eq. (|) with the above, we see that the mapping of the F's versus y 
to the /'s versus ip (or 4'') is area-conserving: the dimensionless scaling variables 
contain a factor kp^, while the scaling functions a factor kp^. In the RFG model 
the /'s display scaling of both the first and second kinds, namely, the display 
superscaling. 



3 Determination of kp and E shift 

In [H, the approach summarized in the previous section was applied to an 
analysis of the usable data involving inclusive electron scattering in the region 
of the quasielastic peak. While medium-energy results were included, the main 
emphasis in that study was placed on the high-energy results from SLAG and 
from TJNAF [^-[^. And, being the first attempt to explore scaling of the 



second kind and superscaling, we chose not to perform an extensive search to 
find the "best" choices of the two parameters involved in the fit, namely, kp and 
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Eshifti but instead selected values that were "reasonable". Now, given the success 
of that previous study and the fact that scaling of the second kind appears to be 
quite well obeyed in the scaling region {ip' < 0), we have stronger motivation to 
produce even better fits and to assess the uncertainties in the fit parameters. 

In the next section we shall place our focus on the transverse scaling function 
fx and it would be desirable to adjust kp and E shift for each nuclear species for 
this quantity. Unfortunately, in the regime where V'' <^ 0, it has not been possible 
to separate the longitudinal and transverse inclusive responses and thus we are 
forced to make our fits to the total /'s. Our hope is that the parameters obtained 
from fits to the total are also appropriate for the individual /l,t- 

The total /'s are very sensitive functions of kp in the region where ip' Q 
and yet it is possible to find values oikp for which the data line up extremely well 
(see, for instance, the high-g/small-?/'' data shown below). The Eshift dependence 
is less critical than is that on kp] nevertheless, by examining the behavior near 
the quasielastic peak it is clear that some shift is needed to move the response 
from the naive peak value oi uj = \Q'^\/2mN (A = r) to where the data require 
it to be. In Table 1 we list the parameters obtained from global fits to the data 

i-in- 



Table 1: Adjusted Parameters 


Nucleus 


kp (MeV/c) 


Esuft (MeV) 


Lithium 


165 


15 


Carbon 


228 


20 


Magnesium 


230 


25 


Aluminum 


236 


18 


Calcium 


241 


28 


Iron 


241 


23 


Nickel 


245 


30 


Tin 


245 


28 


Gold 


245 


25 


Lead 


248 


31 



Clearly the energy shift does not vary too much. Presumably it incorporates 
the separation energy Es-, the mean binding energy of nucleons in the nucleus 
and some global aspects of final-state interactions (for instance, RPA correlations 
which are known to shift the response slightly). While attempts are being made 
to account for the values found here, it is not a simple problem to address for the 
relatively high-energy conditions of most of the current study where relativistic 
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effects are Icnown to be very important, and thus here we hmit ourselves to the 
present phenomenological discussion. 

The values of kp found vary monotonically once "typical" nuclei — say, be- 
yond carbon — are reached. It should be understood that the values given here 
are relative and not absolute: if all values are scaled by a common factor, then 
equally good scaling of the second kind is obtained. The value of 228 MeV/c 
for carbon is typical of other studies and so we have used this to normalize the 
rest. The present fits are done emphasizing the large negative ip' region where 
"contamination" from pion production, 2p-2h MEC effects, resonances and DIS 
are thought to be small and where the basic underlying nuclear spectral function 
is presumably revealed most clearly, in contrast to some previous attempts to de- 
termine kp using the entire response region. Interestingly the values found here 
for heavy nuclei are somewhat smaller than those generally chosen — lead, for 
example, is sometimes assumed to have kp — 265 MeV/c. We beheve that the 
present values are more reliable determinations of the effective kp^s. Note also 
the curious value for lithium, curious because the kp used for ^He is 200 MeV/c 
(see 0). However, this is easily explained if one assumes that ^Li is essentially 
a deuteron (with kp = 55 MeV/c) plus an alpha-particle. Taking the weighted 
mean [{4 x (200)^ + 2 x (55)^}/6]^/^ gives 166 MeV/c, which is very close to the 
fit value of 165 MeV/c. 

To get some feeling for the sensitivity of the fits in Fig. |1| we show the ratio 
Jau/ fc for data from SLAC taken at 6e = 16 degrees and incident electron energy 
Ee = 3.6 GeV — for more discussion of the existing data used 01^lHl see [|l], If. 
The top panel in the figure shows that in the region ip' < the data themselves 
scatter at roughly the 10% level, i.e., scaling of the second kind for these kine- 
matics is satisfied to roughly the 10% level, which is clearly better than scaling of 
the first kind (see [0, Q and figures given below). At positive ip' the ratio moves 
above unity and constitutes the focus of the discussions later in the present work. 
In the middle panel the kp of gold has been increased by 10 MeV/c: clearly the 
fit is very poor in the negative ip' region, indicating that the values of kp given 
in Table 1 are rather finely determined, namely, to only a few MeV/c. In the 
bottom panel in the figure, the energy shift used for gold is increased by 10 MeV 
and again the fit is much worse than the "best-fit" value. So the values of E shift 
given in the table are good to perhaps a few MeV. Finally, note that, in the 
large positive ip' region which will be discussed in depth below, the sensitivity to 
variations in either parameter is much weaker. This stems from the rapid ail-off 
of the responses below = 0, in contrast to the relatively fiat response when 

i)' > 0. 
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4 Transverse Scaling Functions 



From our previous analysis we have seen indications that the longitudinal 
scaling function Jl exhibits superscaling behavior, that is, it not only displays 
scaling behavior of the second kind (as does the total / discussed above), but it 
has scaling behavior of the first kind. Of course, the regime in which this super- 
scaling has been verified is relatively limited, given the difficulty of separating the 
longitudinal and transverse response functions. In practice, using the analysis in 
1 2^ we have made a fit to the combined set of values for the higher momentum 



transfers where scaling of the first kind is seen to occur. The results are shown 
in Fig. g. 

To make any further progress on the problem, it is necessary to make an 
assumption and we do so now: we assume that the so-determined longitudinal 
scaling function shown in Fig. |^ is universal {i.e., superscaling works). Given this 
universal /l we can immediately reconstruct the longitudinal cross section for 
any kinematical condition using the expressions given in Sec. 0, 

= ^hcTMVLGL, (24) 

from this isolate the transverse part of the cross section. 



(25) 



and so obtain the transverse scaling function. 



= V- (26) 

Using this procedure we arrive at the transverse scaling function /y. In Fig. ^ 
we show the results obtained for all kinematics from medium-energy measure- 
ments at 500 MeV and 60 degrees (g = 0.4 GeV/c) to results from both SLAG 
and TJNAF ranging up to g = 4 GeV/c. For ip' < —0.3 we see a reasonable con- 
vergence of the results to a band, although the width of the band is not negligible, 
reflecting (at least) breaking of scaling of the first kind. Since the span of momen- 
tum transfers is so large in the results shown in the figure, we are emphasizing 
the lack of l^*-kind scaling, and focusing on a smaller range of q produces less 
spread, as discussed below. Note also that the region above ip' = —0.3 contains a 
very large spread, that is, a very large degree of scale-breaking. A motivation of 
the present work is to begin to get some insight into the nature of this behavior. 

In Fig. ^ we show only the medium-energy results for fx. Here, at energy 
500 MeV and scattering angle 60 degrees, the momentum transfer varies from 
about 490 MeV/c at ip' = —1 down to about 430 MeV/c at the largest values 
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of ip'- We see on the one hand, that now the band in the negative--?/'' region is 
fairly tight, an indication that the breaking seen in Fig. ^ is indeed mainly due 
to l**-kind scale breaking and not to 2"'^-kind breaking. On the other hand, the 
behavior at positive-'?/'' shows that there one also has breaking of 2"''^-kind scaling 
behavior. Clearly as one proceeds from light nuclei with low kp to heavy nuclei 
with large kp the trend is to increasingly large values of Jt in this region. This 
indicates that the mechanisms that produce the scale-breaking must go as some 
positive power of kp. Indeed, in recent work ^ ^ such breaking of both 1**- 
and 2"'^-kinds due to MEC and correlation effects in the Ip-lh sector has been 
investigated in detail and work is in progress to arrive at a relativistic extension 
of older work p8| in which scale-breaking in the 2p-2h sector was also identified. 

Next, in Fig. ^we show fp for SLAC data at 3.6 GeV and 16 degrees scattering 
angle. At the lowest values of ip' the momentum transfer is roughly 990 MeV/c 
while at '?/'' = 4 it has risen to about 1.7 GeV/c. As the inset on a semi-log scale 
clearly shows, the quality of second-kind scaling behavior at higher g-values is 
excellent in the negative-?/'' region. At positive ip' 2"''^-kind scaling is not perfect, 
although it is only modestly broken. The lower the value of kp the smaller is fp, 
indicating again that the scale-breaking mechanisms go as some positive power 
of kp, as expected. 

Finally, in Fig. ^ we show several sets of data taken at 4 GeV energy and 
various scattering angles. In this single figure one is able to assess scale-breaking 
of both kinds. Each fairly tightly-grouped set of data displays the extent of 
2"'^-kind scale-breaking and clearly the results are very similar to those already 
seen in the previous figure. On the other hand, the various scattering angles yield 
correspondingly different values for the momentum transfer, ranging, for instance 
at ip' = +1 from about 1.2 GeV/c at 15 degrees to 3.9 GeV/c at 74 degrees. For 
each different choice of kinematics a different grouped set of results is obtained, 
indicating again the extent to which scaling of the l**-kind is broken. 

In summary, we deduce from these results for fp that in the region below 
ip' = scaling of the 2"'''-kind is excellent, while scaling of the 1^* is violated, 
although not too badly. In contrast, for ip' > scaling of the 2"'^-kind is also 
broken to some extent, whereas scaling of the l'^*-kind is very badly broken. 
Accordingly, let us next turn to a closer examination of the region above the 
quasielastic peak to see if further insight can be obtained on how the moderate 
scale-breaking of the 2"'^-kind arises. 



5 Ratios of Transverse Scaling Functions 

The results given in the previous section clearly show that several different regions 
are involved when studying the first- and second-kind scaling behaviors of fp. In 
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the region below the quasielastic peak {ip' < 0) at high-g the first-kind scahng 
is reasonably good — this is what is usually called simply y-scaling — and the 
second-kind scaling is excellent. Another example is shown in Fig. containing 
TJNAF data at 4 GeV and 55 degrees (g = 3.4-3.5 GeV/c). Clearly one could 
present the results as a ratio, that is, as [/r('?/'')]2 / where 1 and 2 

denote two different nuclei. For ip' < one would obtain unity with very small 
uncertainties at all but the lowest values of the scaling variable, while for ijj' > 
the ratios deviate from unity. In particular, if 1 light nucleus and 2 ^ heavy 
nucleus, then the ratio rises above unity in the region above the quasielastic peak. 
For instance, in Fig. |^ at ^/^' = +1 the Au/C ratio is about 1.2. Thus, even at 
such high values of the momentum transfer the ratio in the entire region from 
very large negative ip' to very large positive ip' falls within roughly 20% of unity. 

Let us now focus on the region above the quasielastic peak and see if more light 
can be shed on the remaining (modest, i.e., ~20%) second-kind scale breaking. 
To begin with, note that the scaling functions /, //, and /t have been constructed 
with the quasielastic responses in mind: the denominator in Eq. ([13D contains the 
elastic eN form factors weighted by Z and N, the proton and neutron numbers, 
respectively. In the region below the quasielastic peak we expect the dominant 
contributions to the response functions to arise from the tails of the QE response 
and accordingly removing such a weighted elastic eN cross section makes sense. 
However, in the region above the quasielastic peak we do not expect things to 
be so simple. Certainly resonance excitation and pion production play an impor- 
tant role and, at asymptotically high values of momentum transfer, one expects 
DIS to take over. Because of this one expects the longitudinal and transverse 
responses to have quite different characters — for instance, the first important 
new contribution beyond QE scattering in entering the ip' > region is expected 
to come from N —>■ A excitations on nucleons moving in the nucleus, and this is 
(1) essentially transverse and (2) isovector. For this reason we have extracted the 
transverse response using the superscaling assumption for the longitudinal con- 
tributions. Then, focusing on the transverse results, we can attempt variations 
on the theme that motivated y-scaling in the first place. Specifically, instead of 
weighting protons and neutrons with G\jp and G\.j^, respectively, we can make 
a scaling function that varies only with total mass number A = Z + N and not 
Z and individually — this is more in the spirit of what is done for the ratio 
considered when studying the EMC effect. For instance, if an isovector excitation 
such as A^ — A dominates in the response, then equal weighting of protons and 
neutrons should be better that the weighting provided using the elastic magnetic 
form factors of the nucleons. Accordingly let us consider 

- _ kpTjT 
^~ AaMVT{l) Gl' 

12 



(27) 



where Gd is the dipole form factor (a compromise between the dependences of 
Gmp and Gmu — its form is not important here as it will cancel in constructing 
the ratio defined below). This provides an alternative for Jt defined above in 
Eqs. ([T5|j2^) . The two differ roughly by a factor 



ZGljp + NG\j^ 
AGl 

^ 2^^ 2 

5.73 [1 



0.36(iV - Z)/A] 



(2J 



which is a relatively slowly- varying function of (Z, N). Such an approximation 
will allow us to interconnect the quasielastic and DIS analyses in an approximate 
way (see below). 

With these definitions, let us next consider two different nuclei have charge 
and neutron numbers (Zj, Nj) with ? = 1, 2. Labelling as above with nuclear 
species numbers 1 and 2, we consider the ratio 



Pl2 



(29) 



With -fil^ 



Aii^//^2^ we then have that 



Pl2, 



(30) 



where the last approximation holds to the degree that 

_2 ^ 1 - 0.36(iV, - Zi)Mi 
1 _ o.36(Ar2 - Z2)M2 



(31) 



is nearly unity across the periodic table. For instance, in going from a nucleus 
such as carbon {Ni = Zi) to gold (A^'2 = 118, Z2 = 79) one has 71^2 — 1-08. Thus 
the ratio of the /'s will be about 8% larger than the quantity pi2 for gold/carbon 
and even closer to unity for ratios involving nuclei closer together in the periodic 
table. 

From this simple analysis we can see that some of the rise discussed above 
for the ratio [/r(V'')]2 / [/t(V'')]i ^^^e region ip' > occurs because of the differ- 
ent weighting of protons and neutrons in regions dominated by processes other 
than quasielastic scattering. Two examples of the ratio pi2 are shown in Fig. ^. 
These have been obtained by extracting the transverse contributions from the 
experimental cross sections for each nucleus using the superscaling assumption to 
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remove the longitudinal contributions, as discussed in the previous section, calcu- 



lating frii^') for each nucleus using Eq. (^), fitting the results and finally using 
Eq. (^) for specific pairs of nuclei. The upper panel in the figure corresponds 
to momentum transfers of about 1.2 GeV/c (near ip' = 1) to 1.7 GeV/c (near 
ip' = 4), whereas the lower panel ranges from about 2.1 GeV/c (near ip' = 1) 
to 2.4 GeV/c (near ip' = 2.5). Thus, as the momentum transfer grows the ratio 
seems to stabilize at roughly 5-10% above unity (no error bands are shown in 
the figure, but typically the uncertainties are ±2%). The analysis in terms of pu 
accounts for some of the 2"''^-kind scale breaking in the region above the quasielas- 
tic peak; however, it still does not explain the full effect, but only perhaps half 
of the 20% discussed above. To get more insight we need to invoke a model of 
the transverse response in the region ip' > 0. We proceed in the following two 
subsections by examining the resonance region and DIS separately. 



5.1 In the Resonance Region 

Let us begin by focusing on the region in which we expect the excitation of 
resonances and meson production to become relevant, in addition to the tail of 
the quasielastic response. In Fig. || this corresponds to ip' ~ 1-2 (upper panel) 
and somewhat smaller (lower panel); see below where the scaling variable ip^ is 
introduced for a more quantitative measure of where a given baryon resonance 
has its peak. In |^ an extension of the RFG model was discussed. Instead of 
building the model from elastic scatterings of electrons from nucleons in this work 
the inelastic N A transition was considered and again a relativistic Fermi gas 
model constructed. Of course the ideas can be extended even further to inelastic 



excitation of any baryon resonance, that is, with any mass m^,. Following p9 



the RFG hadronic tensor for the production of a stable resonance of this mass is 
given by 



W, 



N 



kp 



dp 



E{p)E4p + q 
x6[E,{p + q)~E{p 



0{kF-p)ff,^{p,p + q 



UJ 



(32) 



where N' is the number of protons or neutrons (the total nuclear response should 
be Z times the result for target protons plus times the result for target neu- 
trons). As usual E{p) = Jp^ + mj^ and now also E^{p + q) = J{p + qY + ml. 



Here fav is the single-nucleon inelastic hadronic tensor: 



Uu{p,P + q) 



-Wi[T 



9iMu 
1 



P -^Q 



p„ 



g2 



(33) 
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where wi^2 are the analogs of Wi^2 for elastic scattering — the w's for the ^ A 



case are discussed in [gg] ; however, for the present purposes it is sufficient simply 
to know that they are functions only of the 4-momentum transfer r. 

For the strict RFG model one assumes that an on-shell nucleon moving in 
the initial-state momentum distribution is struck by the virtual photon and an 
on-shell baryon resonance of mass m^, is produced. Of course this is a very 
over-simplified model (the struck nucleons are not on-shell, the single-nucleon 
process is not all, the final state is not a stable baryon but a decaying resonant 
or non- resonant state, etc.), although it does provide some insight into a possible 
scaling violation mechanism. In particular it can shed some light on when such 
scale breakings do or do not appear to be very large; that is, it helps to define 
characteristic kinematic conditions where the QE-type scaling could be expected 
to evolve into another type of scaling behavior. 

Specifically, using the N ^ A RFG model of P^] and forming the ratio pi2 



defined above one typically obtains results above unity in the ip' region lying 
between the quasielastic and "A" peaks. This is to be compared with the experi- 
mental ratio, shown in Fig. ^ For the upper panel in the figure the — > A RFG 
peak occurs at ip' = 1.3 for carbon (1.2 for gold), while for the lower panel the 
peak occurs at ip' = 0.72 for carbon (0.67 for gold). The model yields a typical 
ratio Pau/c in the region between the QE and — A peaks of approximately 
1.08. For the conditions of the lower panel the model yields a ratio that is closer 
to unity, approximately within 4-5% of unity in the region between the QE and 
iV — i> A peaks. Thus we see, even with this over-simplified model, that the re- 
maining few percent deviation from unity of the ratio pi2 is reasonably compatible 
with resonance production in the z/'' > region corresponding to excitation of 
the A. 

Some further insight can be gained about why the ratio behaves as it does by 



examining some of the specifics of the modeling. In fact, in a scaling variable 
emerged naturally from consideration of the RFG model for the "A" region. Here 
we generalize this to the scaling variable that occurs in this model for inelastic 
single-nucleon transitions (m^r — ^ m^). With 

P* = l + -, (34) 
r 

where = {ml — m%)/4:m% > 0, we have the analog of Eq. (^: 

^* = , (35) 



liF^(l + Xp.)T + K^T{Tpl + l) 

Futhermore, following the discussions in Sec. |^, one may shift the energy to obtain 
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the analog of Eq. 



^* = ^ I (36) 



^(1 + Ayjr' + K^T'ir'p',^ + 1 

where p'^ = 1 + 13^/t' . If = mTv, then = 0, which imphes that = 1 
and the above generahzed scahng variable yields the standard results in Eqs. (^ 
and (^. The RFG model for inelastic excitation of a stable resonance of mass 
m^, has scaling behavior of both the first and second kinds in the variable tp^. In 
particular, if a universal transverse scaling function Jt underlies both quasielas- 
tic scattering and the excitation of some given resonance from nucleons in the 
nucleus, then we expect the former to occur in the total transverse cross section 
as fxii^') weighted with the usual eN elastic cross section and the latter as /t(V'I) 
weighted by the N— >-N* inelastic cross section. Thus we expect (at least) that the 
scaling behavior involves the way /t(V'') varies with tp' — it scales — versus the 
way /t('?/'*) varies with ip' — note: not how it varies with ip^ in which it scales. 

Relationships can be established between '^|J^, and ip (or ipj and ^p' using the 
primed kinematic variables as above — in the following for simplicity we consider 
only the unshifted variables): 



(37) 



where 27^ = [(1 + r)(l + p^r)]^/^ ^ (1 + P*t). Substituting for p^ one can write 

-1/2 



7* 



2r(l + r) 



2Jt(T+V) 




(38) 
(39) 



where the limiting behavior is for r ^ 1 for fixed {i.e., fixed m^,). Under these 
conditions 7* falls as [r(l + t)]~^/^. Thus, if one keeps m^, fixed, focusing on a 
particular resonance, for instance, and examines the behavior of ip* as a function 
of T and ip it is clear that as r becomes very large 7* ^ and therefore ip^, ^ ip, 
so that there will be no distinction between the two scaling variables. Indeed, to 
get an idea of how large r must be for this to occur, suppose that in addition one 
has \riF'ip\ 1) then a simple approximation for %p^ is the following: 

^ V (40) 

r]F\r{l + T) 
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So the value of r that characterizes the coalescence of ip* and is tq = ([1 + 
(2/5*/?7f)^]^''^ — l)/2. For instance, if m^, = and /c^ = 230 MeV/c, then one 
has To = 0.4, corresponding to \Qq\ = 1.4 (GeV/c)^. As the mass of the resonance 
grows then so does the characteristic tq, and the coalescence is postponed until 
larger values of r. 

Given these relationships between the various scaling variables, let us now 
return to the analysis of 2"'°'-kind scale breaking in the resonance region. As noted 
above we expect some linear combination of /t(V^') and /t(V'*) with weightings 
according to the elastic and inelastic eN cross sections to be involved in the region 
between the quasielastic peak {ip' = 0) and the peak of the specific resonance of 
interest {ipl = 0). As a function of ip' the quasielastic contribution is assumed 
to exhibit 2"''^-kind scaling, whereas the N— i>N* contribution does not, because of 
the dependence on kp in the expressions above. Consider for simplicity the result 
in Eq. (^D|) . With r fixed and considering a given value of and hence of (3^, we 



see that ip* is less than ip by the offset P*/[i]f\/t * (1 + r)]; so is some positive 
number whereas (being between the QE peak and the resonance peak) ip^, is some 
negative number. For a small kp the offset is large, whereas for a large kp it is 
small. That is, for kp small ip^ is further away from zero than when kp is large. 
Or, said a different way, the difference 



is positive if kpi > kp2, and so ip^{kpi) is closer to zero in this case than is 
V'*(^F2)- Since the inelastic response peaks when ip,, = 0, this implies that in the 
region between the QE and N^N* peaks the case with the larger kp should be 
larger than for the lower-kp case. Indeed this is just what we observe from the 
data and, in fact, if the offset given above is computed for typical conditions, the 
degree to which scaling of the 2"'^-kind is broken is found to be quite compati- 
ble with those results. These ideas have natural extensions to a distribution of 
resonances between each of which the same type of phenomenon will occur. The 
actual magnitude of the effect of course depends on the cross section weighting 
of each contribution and thus the present arguments should only be regarded as 
qualitative ones — a more quantitative approach is presently being pursued. 

In summary, in effect, the kinematic offset from the QE scaling variable which 
occurs in the scaling variable that enters naturally when considering baryon ex- 
citation from nucleons in the nucleus alone appears to be capable of explaining 
most of the breaking of 2"'^-kind scaling in the resonance region. 




(41) 
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5.2 In the DIS Region 



Turning now to the DIS region, let us recall the conventional analysis of x-scaling 
(see, for example, ||3^). Noting that the quantity z/ employed in DIS analyses is 
called uj in our work, as is usual, and employing the above dimensionless kinematic 
variables we can write 

x=^- (42) 



and see that when < we have x > 1 and the reverse; when = we have 
X = 1. 

Using the fact that x = rjX one can construct expressions that directly relate 
X to the scaling variable if) and Fermi momentum (via ^p) — in particular one 
useful form is the following: 



X 



2 




(43) 



Or alternatively, if one wishes to work with the Nachtman variable ^ (not to 
be confused with used in this and other RFG-motivated studies), then again 
various useful relationships exist. In particular. 



2(k- A) 



H 2 



(44) 
(45) 



Moreover, it is possible to incorporate the energy shift E shift as above by replacing 
all quantities on the right-hand sides of Eqs. (|^3|-^) by their primed counterparts 
to obtain shifted variables x' and 

Let us examine the x-variable a little more closely. As written in Eq. ( ^3]) 
it is a function of r, (and hence kp) and if). The last two enter only in the 
combination z e 



ip/'lil) = ripijj/2, since 



1 

X 



Vl + z^ + z 



+ 



1 



1 + - - 1 } 2zVl + 



(46) 



Thus, as T 



oo one finds that 

1 1 



Vl + z^ + z . (47) 

X X'^ ^ J 

This behavior holds for r ^ Tcru, where, by comparing x with x°° in the regime 
where ip > 0, one finds that 



Tcrit 



Z\/\ + Z^ 

VTTz^ + z 



< 1/4. 



(4J 
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That is, for ^ 1 (GeV/c)^ one has that x is well approximated by x°°. 
In this regime x becomes a function only of z and thus DIS scaling in x {i.e., 
independence of r) and l'^*-kind scaling in ip will occur together. Furthermore, 
the above expression for yields 



= 1 + 2^ for < V2/r]F 
^ Az^ for ^ > V2/r/i.. 



X' 



(49) 
(50) 



The "asymptotic regime" then occurs when \Q'^\ ^ 1 (GeV/c)^ and x is small 
enough that xjj ^ V^/r]F (typically = 5-6), and under these conditions 



X 



(51) 



Thus in the asymptotic regime x oc 1/kp. 

Re-expressing what we have seen in the analyses so far: throughout the region 
X > 1, near x = 1 (essentially the QE peak) and down to about x = 0.3-0.4 
(corresponding to ip' = 3-4 in Fig. H) we find that the ratio pi2 is close to unity 
and that its deviation can be reasonably accounted for by the modest second- 
kind scale-breaking expected for resonance excitation. Now we wish to explore 
the higher- 1 I , small-x region in which DIS is expected to be the dominant 
reaction mechanism. 

We begin by providing some "translations" between the nomenclature and 
conventions employed for most studies done within the context of nuclear physics 
and those done in particle physics studies. Specifically, what are usually called 

and ax in the latter are given in terms of the responses Wi^2 and hence via 
Eqs. (|12|) of Rl,t by the following (here K is an arbitrary factor conventionally 
taken to be q or sometimes u — uj{ip = 0)): 



K 
K 



Wi 



Their ratio R is then 



—W2 - Wi 

T 



T 



not to be confused with a different ratio 



Rt 2t 



(52) 
(53) 

(54) 
(55) 
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As in the previous discussions, we wish, to the degree that we can, to isolate 
the transverse part of the inclusive cross section. In the DIS regime we must 
depend on the limited information available for the ratio R defined in Eq. 
STI, |32i- We then construct 

a a , ^ 

(56) 



1 + 



1 + 



and from this the transverse scaling function J't and the ratio pi2, as above. In 
addition to the results shown in Fig. ^ which are not really in the DIS region, 
we show in Fig. |^ some of the (limited) information available at higher- — 
apparently no information on individual eA cross sections is available at very 
high momentum transfers and only the usual EMC ratio at constant x has been 
measured. 

Scaling in the DIS regime is usually discussed in terms of the quantities Fi = 
ttinWi and F2 = VW2. The so-called EMC ratio is chosen to be [F2/ A]2/[F2/ A]i, 
where, as above, 1 and 2 label a given pair of nuclei. When this is plotted versus 
X at very high \Q'^ \ it is seen to be independent of and to exhibit the so-called 
EMC effect. Now our previous analyses and the above discussions show that if 
one chooses for the abscissa on a plot to use the scaling variable ip' , scaled by the 
Fermi momentum, then one should also incorporate the charge of variables into 
the ordinate in the plot by multiplying hj kp. Explicitly this was the procedure 
used with ip' = y/kp (see Eqs. (|,^) and f = kp x F (see Eqs. (p2|,p3|)). Since 
here we always display results as functions of tp', we are then motivated also to 
employ a modified "EMC-type" ratio involving the structure functions for the 
two nuclear species 1 and 2, but now with factors of kp for each: 

nis _ [kFF2/A]^ 
[kpF2/A\^ 

Let us now provide the connections between the QE-type ratio and the EMC- 
type ratio. These may be related via 

pf/^ = X12P12, (58) 

where the factor X12 in turn may be decomposed into three factors, 

X12 = R12G12K12. (59) 

The first of these, 

arises simply because different parts of the total cross section are conventionally 
chosen in forming the two ratios. When comparing two different nuclear species 
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at constant ip' it must be remembered that x for species 1 is not the same as x 
for species 2, but that the two x-values are shghtly different. In fact, R12 is quite 
close to unity for the range of kinematics under study in this work. The second 
factor, 



G 



12 



[G 



D\2 



(61) 



stems from the fact that the QE ratio is defined in Eq. (|29|) via Eq. (|27[) to have 
the typical dipole form factor behavior expected for the transverse part of the 
elastic eN cross section, whereas the DIS ratio is not — in the latter we expect 
e-quark scattering to begin to dominate. Again, the fact that we display results 
versus ip' for experiments where incident electron energy E^, and scattering angle 
6e are fixed means that for the two nuclear species the 4-momentum transfers 
are slightly different. Indeed, G12 is typically quite close to unity for the range 
of kinematics under study. Note that, if sufficient information were available to 
allow one to specify \Q'^\, rather than have it vary as here (excepting the case 
shown in Fig. P), then this factor would be irrelevant. 

The important factor for our purposes is the third one. 



K 



12 



(62) 



which contains all of the kinematic behavior involved in relating the two types of 
ratios. At the quasielastic peak this is unity (for simplicity Eghift is neglected in 
these arguments), whereas above the peak it in general differs from unity and so 
must typically be taken into account. To get some feeling for how it enters in the 
extreme DIS regime where r ^ 00 and x ^ 1, note that there one has 



-"■12 



'x^ 










2 


'x^ 










1 



(63) 



Thus in extreme-DIS kinematics the ratio defined in the present work in Eq. ( p9D 
for use in QE scattering differs from the conventional EMC-type ratio essentially 



by the factor [a; ]2/k ]i in Eq. 



From Eq. (^) above we see that this 



immediately implies that in the asymptotic regime 



^12 = 



kp 



(64) 



Returning to Fig. ^, we see from the lower panel that scaling of the second kind 
in this kinematic region is not as good as it was for the resonance region and below. 
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Multiplying pu by the ratio of the x's as in the upper panel in the figure yields 
results much closer to unity, as expected. Namely, the upper panel corresponds 
to the EMC behavior when results are displayed versus ip'. Indeed, uncertainties 
are not shown in the figure as they are hard to quantify; what is clear is that the 
results in the upper panel are essentially consistent with only small deviations 
from unity. Finally, note that from the above arguments asymptotically we would 
expect factors of [kp]^ / [kp]2, which from Table 1 are 1.15, 1.25 and 1.33 for Al/C, 
Fe/C and Au/C, respectively. While not unreasonable in fact these are somewhat 
further from unity than the actual ratios in the lower panel in the figure. Indeed, 
while \Q^\ is reasonably large here (r = 1.42 > Tcru in Eq. P8|)), x is not very 
small, ranging as it does from roughly 0.7 at the left-hand side of the figure 
down to approximately only 0.25 at the right-hand side. Or, in other words, the 
asymptotic regime where ip' is sufficiently large for Eq. (|50|) to be valid is only 
barely reached at the right-hand side of Fig. ^ and the behavior seen here refiects 
the fact that these data lie somewhere between the quasielastic and asymptotic 
regimes. 

6 Discussion and Conclusions 

In undertaking this study we have begun by refining our previous analysis of 
the kinematic region lying below the quasielastic peak. We had seen in our past 
work that scaling of the second kind is excellent in this region and that values of 
kp and Eg^ift can be determined for each nuclear species for which medium- to 
high-energy data exist. In the present work we have taken these determinations 
a step further and evaluated how sensitive the scaling function is to variations 
around the best-fit values oi kp and E shift-, finding that the former is very well 
determined (relative to an overall multiplicative factor), whereas the latter is less 
important, as found in our previous studies. 

We have then employed the hypothesis of longitudinal superscaling to extract 
the longitudinal contributions from the total inclusive response and so obtain the 
transverse response and scaling function. It should be stressed that the kinematic 
regime in which the superscaling of /l is tested experimentally is unfortunately 
rather limited and hence the extraction of /t over a wide range of kinematics is 
contingent on having fi superscale. However, it is important to note that espe- 
cially at high g, where the L/T ratio becomes small, the cross section is known 
to be dominated by the transverse response. Our focus in the present work is 
this region and thus the extraction procedure should be reasonably good here. 
Proceeding with this as a basic assumption using the fi determined by what TL- 
separated data can be relied upon, we deduce fx for an extended range of kine- 
matics. We expect that breaking of scaling of both the first and second kinds will 
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be more prominent in fx than in due to the well-known reaction mechanisms 
— for instance, pion production, N ^ A excitations and 2-body meson-exchange 
current effects, all of which are expected to be transverse-dominated — and ac- 
cordingly have placed our focus on Jt- This is not to say that /x, cannot also 
have interesting scaling-breaking behavior — for instance, final-state interaction 
effects such as those arising in the random-phase approximation can differ in lon- 
gitudinal and transverse, isoscalar and isovector channels — however, isolating 
/l except as we do here via the universal superscaling approach has proven to be 
very difficult and thus the focus on /t is inevitable at present. 

We have selected two regimes for the present extended study, (1) the reso- 
nance region above the quasielastic peak where we expect N ^ A and N* 
contributions to be important, and (2) still higher inelasticity where DIS takes 
over. Unfortunately, there are very few cases for the latter in which eA cross 
sections are available for a range of nuclei; most data involve ratios involving a 
pair of nuclei for x held constant, not what we want, namely, the individual cross 
sections as functions of \Q'^\ and scaling variable ip'. 

What we observe is that scaling of the first kind is clearly violated when 
one proceeds above the QE peak (■0' > x < 1), whereas scaling of the 
second kind is better. Indeed, for sufficiently high momentum transfers it is 
excellent for ip' < ^ x > 1 and typically is broken only by 10-15% up through 
the resonance region. Moreover, from simple modeling and even from simple 
kinematic arguments involving the introduction of the scaling variable that 
occurs naturally when discussing the electroexcitation of resonances, this modest 
second-kind scale-breaking is reasonably accounted for, with very little further 
scaling violation left to explain. This does not leave very much room for the 
other reaction mechanisms that are known to violate scaling of the second kind. 
For instance, we know that at least some of the meson-exchange effects depend 
on kp very differently than does the usual 1-body QE contribution (typically in 
a way that breaks scaling of the second kind by contributions proportional to 
kp), and thus the absence of significant contributions with this behavior limits 
how much of a role they can play. In only a very limited number of cases has 
it been possible to carry through modeling of MEC effects with the requisite 
relativistic content and at present a concerted effort is being made to extend 
past non-relativistic treatments of MEC and interaction effects to the kinematic 
regime of interest and thereby to test these ideas in more depth. 

Finally, in the DIS region we see that there is clearly more breaking of scal- 
ing of the second kind, indicating (as expected) that the reaction mechanism 
is becoming different, presumably e-quark physics rather than eN physics. The 
cross-over from the resonance region, where second-kind scale-breaking is reason- 
ably small and can be explained using simple arguments, to the DIS region, where 
such does not appear to be the case, may in fact be a relatively clear indicator of 
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the shift from hadronic to QCD degrees of freedom. 
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Figure 1: The ratio of Jau over fc for energy 3.6 GeV and scattering angle 16 
degrees. The top panel shows the ratio using [/cf]au ~ 245 MeV/c, [EghiftjAu = 25 
MeV and [kF]c = 228 MeV/c, [EsMftlc = 20 MeV, the "best fit" values. In the 
middle panel [/cfIam has been increased by 10 MeV/c, while in the bottom panel 
[EghiftlAu has been increased by 10 MeV. 
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Figure 3: (Color) Transverse scaling function Jt for a wide range of nuclei and for 
kinematics ranging from medium energies (500 MeV and 60 degrees) to high energies 
(up to 4.045 GeV and 74 degrees). The longitudinal response has been removed using 
the superscaling assumption discussed in the text. The fit parameters are listed in 
Table 1. 
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Figure 4: (Color) As in Fig. ^, but only for medium energy data (500 MeV and 60 
degrees) . 



31 




Figure 5: (Color) As in Fig. but only for data at 3.6 GeV and 16 degrees. 
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Figure 6: (Color) As in Fig. but for data at 4.045 GeV and angles ranging from 
15 to 74 degrees. 
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Figure 7: (Color) As in Fig. but only for data at 4.045 GeV and 55 degrees. 
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Figure 8: The ratio pi2 defined in the text for two kinematic conditions, 3.6 GeV 
with 16 and 30 degrees. 
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Figure 9: The ratio pi2 defined in the text for = 5 (GeV/c)^ and the same 
ratio multiphed by (a:2)^/(xi)^ for pairs of nuclei 1 and 2. 
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